THE GRILLAKIS-MACHEDON-MARGETIS SECOND ORDER 
CORRECTIONS TO MEAN FIELD EVOLUTION FOR WEAKLY 
INTERACTING BOSONS IN THE CASE OF 3-BODY 
INTERACTIONS 
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Abstract. This note shows that, with a little modification, the results in [l] 
hold in the 3-body interactions case. 



1. Introduction 

Recently, a second-order correction to the usual tensor product (mean-field) 
approximation for the Hamiltonian evolution of a many-particle system with 2- 
body potential was found in Grillakis-Machedon-Margetis|l . Here is their main 
result, 

Theorem 1. [T] Let (f> be a smooth solution of the Hartree equation 

(1.1) ^ + A«£ + (</*|0| 2 )^ = O 

at 

with initial conditions 4> and potential v' being even. Assuming the following: 

(1) Let k(t,x,y) E L 2 s {dxdy) for a.e. t, solve the Grillakis-Machedon-Margetis 
equation 

(1.2) (iu t + ug T + gu - (I + p)m) — (ip t + [g,p] + um)(l + p)~ 1 u , 
where products in equation mean compositions of operators, and 

u(t, x, y) := sinh(fc) := k + + • ■ ■ , 

5(x — y) + p(t, x, y) := cosh(fc) := S(x — y) + Tjy^ + • • • > 

g(t,x,y) := -A x S(x - y) - v(x - y)(j){t,x)4>{t,y) - (v * |0| 2 )(i, x)S(x - y) , 

m(t, x, y) := v(x - y)(f>(t, x)<f>(t, y) . 

(2) Functions f(t) := \\e B [A,V']e- B VL\\ r and g(t) := \\e B V'e- B n\\yr are locally 
integrable, where 

(1.3) B(t) = i J (k(t, x, y)a x a y ~ k(t, x, y)a* x a* y ) dxdy 
and T the Fock space defined in [1] . 



The author thanks Matei Machedon for discussions related to this work. 
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(3) J d(t, x, x) dx is locally integrable in time, where 
(1.4) d(t,x,y) = (z sinh(/c) t + sinh(fc)g T + gsmh(k)) sinh(fc) 

— (i cosh(fc) t + [g, cosh(fc)]) cosh(fc) 



— sinh(fc)mcosh(fc) — cosh(fc)msinh(fc) . 
Then, there exist real functions \o, Xi such that 

(1.5) \\ e -VN A(t) e -B(t) e -t J*(Nxo(s)+Xi(s))ds n _ e ttH' N £ - VNA(0) qii ^ 



< Iof( s ) ds , Io9(s)ds 



(N N 
where f2 = (1, 0, 0, • • • ) € J 7 is the vacuum state, 
(1.6) A(<£) = a(0) - a*(<£) 

with a and a* are the annihilation operator and creation operator, and the 
Hamiltonian 

H' N = J a* Aa x dx + J v' {x — y)a* x a* y a x a y dxdy 

To prove Theorem [T] the key is to notice that, for <p and k satisfying the assump- 
tions, there is some "nice" Hermitian L such that 

and 

\\LSl\\r < N-^ 2 \\e B [A, V\e- B Ciy + N- l \\e B Ve- B n\\ F , 
where 

Then applying the energy estimate to 



would give the estimate 11.51 

With the same set up, this result can be extended to the case of 3-body interac- 
tion potentials. 

2. Main Result 

We consider the Fock space Hamiltonian Hn ■ T — > T with the 3-body interac- 
tion potential defined by 

(2.1) H N = J a* x Aa x dx + -^-^ J v(x - y,y - z)a* x a*a* z a x a y a z dx dy 
in correspondence to the Hamiltonian 



n 1 n 



H N . N - ^ , N2 

j=l i<j<k 
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We assume that v is symmetric in x, y, and z. As shown in (1.2) and (1.3) in 
Chen-Pavlovic[2j, every translation invariant 3-body potential can be written in 
the form v(x — y,y — z). However, I don't assume v ^ as in |5]. 

Theorem 2. With the same assumptions as in theorem]]] with H' N replaced by Hn, 
V by V in formula \2.1\ We change equation \1.1\ to 



d , . , 1 

i 

and 



(2.2) i^> + + ^ / v(x - y, y - z) \^(y)\ 2 \cf>{z)\ 2 dydz = 



g(t,x,y) := -A8 l2 - ( J v 12 , 23 |</> 3 | 2 dzjfafa Vl2 > 23 \^ d V dz j s i2 

m(t,x,y) := (J ^12,23 |^ 3 | 2 ^^0i^2 

where ^12,23 l^l 2 is an abbreviation for the product v(x — y,y — z) \(j>(z)\ 2 , 
Further, we assume that the functions 

h(t) :=\\e B [A, [A,V]]e- B ny 

and 

i(t) :=\\e B [A, [A, [A,V]]]e- B ny 

are also locally integrable. 

Then, there exist real functions Xo> Xi such that 

^ e -VNA(t) e -B(t) e -i f*(N Xo (s)+Xi(s))ds n _ e itH N e -^fNA(0) n 

< Jo fi s )ds Jl g{s)ds /J h(s)ds J* i(s)ds 
~ 6iVi 67V 2 127V 367V^ 

Proof. Similar to the proof of Theorem [TJ □ 



Remark 1. Recall that, the equation for k is \1.2\ 

(2.3) (iu t + ug T + gu - (I + p)m) = (ip t + [g,p] + um)(l +p)~ 1 i 
A derivation of this equation will be given in Section^ 

Remark 2. In Section^ we will show that the L in this case is 

(2.4) L =H Q + J Vi2,23 I ^3 1 2 4> 2 4> 1 a x a y dxdydz+ 

fi2,23 \4> 2 \ 2 |0 3 | 2 a* x a x dxdydz 
1 



Bjr-B 



d(t,x,y)a* y a x dx + —^e Ve 
' Te B [A ,V]e- B + -±-e B [A,[A,V}]e- 



6JV 3 / 2 1 ' J 12JV 
' * B [A,[A,[A,V}]}e- B 



36N2 

where d(t,x,y) is the same as equation] 1. 4 
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3. The New Hartree equation 



In this section we derive the Hartree equation 12.21 for the one-particle wave 
function <p as needed in Theorem [5J 

Lemma 1. The following relations hold, where A denotes A(cj)), and A , V are 
defined by formulas \1.6\ and \2.1\) : 



3 J v(x — y,y — z)((f>(x)a* y a z a x a y a z + 4>(x)a x a*a* z a y a z )dxdydz 
[A, [A, V]} 

6 J v{x - y,y - z)((j>{x)(t>{y)a* z a x a v a z + 2(f>(x)(f>(y)a* x a* z a y a z 

+(j)(x)(j)(y)a x a y a* z a y a z )dxdydz 

+6 / v(x — y,y — z) \<fi(x)\ 2 a y a*a y a z dxdydz 



[A, [A, [A,V]}] 

= 36 J v(x - y, y - z) \4>{x)\ 2 (4>(y)a* z a y a z + 4>(y)a*a* z a z )dxdydz 

+6 J v(x -y,y - z){4>{x)4>(y)(j){z)a x a y a z + (j)(x)<f)(y)<j>(z)a x a*a*)dxdydz 

+18 J v(x - y,y - z){(j){x)(j){y)(j){z)a* z a x a y + (f>(x)(f>(y)4>(z)a x a*a z )dxdydz 

[A, [A, [A, [A,V]}}] 
= 72 J v(x - y, y - z) \(f>(x)\ 2 {(j){y)<j){z)a y a z + <f){y)cj){z)a* y a* z )dxdydz 

+ 144 J v(x — y,y — z) \(f)(x)\ 2 (f)(y)(f>(z)a z a y dxdydz 

+72 v(x - y,y - z) \(j>(x)\ 2 \4>{y)\ 2 a* z a z dxdydz 



[A, [A, [A, [A, [A, V]]]]] 
= 360 J v{x-y,y- z)\cf)(x)\ 2 \<j)(y)\ 2 {(t){z)a z + (f){z)a* z )dxdydz 

[A, [A, [A, [A, [A, [A,V]}}}}} 
= 720 / v(x - y,y - z)\cj)(x)\ 2 \(j)(y)\ 2 \4>(z)\ 2 dxdydz 



Proof. This is a direct calculation. □ 
Now, we write * x (i) = e VNA(t) e itH e -VNA(o)Q for we carry out tne ca i_ 

culation in the spirit of equation (3.7) in Rodnianski-Schlein [3]. 

Proposition 1. Say <fi satisfies the Hartree equation 

d If 2 2 

i Ql ( f ) + A( l )+ 2^ v '" x ~ v,y ~ z ^^ y ^ d y dz = 
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then ^i{t) satisfies 

= (h n + ~[A, [A, [A, [A,V]}}\ + l -N-^[A,V] 



+ ^N-^A, [A, V}} + ij\T * [A, [A, [A, V}}} 



N 
~3 



/ v(x-y,y-z)\<t>(x)\ 2 \<j>(y)\ 2 \<f>(z)\ 2 dxdydz)il> 1 (t) 



Proof. Applying the formulas 

(e- C{t) )=C+^[C,C} + ^[C, [C,C}] + 



dt 



e c He- c = H+[C,H} + ^ [C, [C, H\] 



2! 

to C — vN A we have 



(3.1) ~Vi(t)=-Mi, 

1 at 

where 

Ll = -( ILe^w) e -^Mt) +e VNA(t) HNe -VNA(t) 



i \dt 

= \ [n^A + ^[A,A]^+H N + NV*[A, H q ] + ^ [A, [A, H ]] + \ (n~V 2 [A, V] 
+ ^[A, [A,V]] + ^-[A, [A, [A, V}}} 
+ [A, [A, [A, V]]]] + ^-[A, [A, [A, [A [A, V]]]]] 

+ ^[A,[A,[A[A[A[AV}]}}}] \ . 

The Hartree eauatior l2.2l is equivalent to setting 

(3.2) -A + [A, #0] + ~ ^ [A, [A, [A, [A, [A, V]]]]] = . 

1 6 5! 



Or more explicitly, equation 13.21 is 

j / 2 2 ^ / 2 2 

a(icj) t +A(t)+-(j) / vi2,23 |0 2 l 1031 dydz)+a*{i(t) t +A(j)+-(f> / ^12,23 |0 2 l 1 03 1 dydz) = 

because of Lemma [T] and the fact that [A x a x , a*] — (A5)(x — y). 
Thus 

- [A, A] + [A, [A, H }} + I i [A, [A, [A, [A [A, [A V}}}}}} = , 
« 5! 
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i.e 13. II simplifies to 

(3.3) = fa + [A, [A, [A, V]]]] + l -N-^[A,V] 

+ ^N-^A, [A, V}} + [A, [A, [A, V}}} 



N 
~3 



^ I v(x-y,y-z) \0(x)\ 2 \<j>(y)\ 2 \<j>(z)\ 2 dxdydz ) * x (t) 



□ 



Write the last term as 

N 
~~3 



^- / v{x-y,y- z) \<j){x)\ 2 \4>{y)\ 2 \(j>(z)\ 2 dxdydz := -iVx > 



then the first two terms on the right-hand side of 13. 31 are the main ones we need to 

consider, since the next four terms are O (-fis), O f^s"J > O (^) ,and O (^r)- 

In order to kill the terms involving "only creation operators" i.e aja* in jrg^, [A, [A, [A, V]]]], 
we introduce B fsee ll.3[) and let 



as in pp. So we have 

1 a 



where 

1 / a 



i 



at 



L= ll -e B )e- B + e B L ie - B 



= L Q + ^~ 2 e B Ve- B + \N-^e B \A, V] e - B + ^e B [A, [A, V]] e - B 
+ ^N-h B [A, [A, [A, V]]]e- B - N Xo , 



and 



(3-4) L Q = 1 f^-A e- B + e B (h + ~[A, {A, {A, {A, V}}}} ) , }; 



4. Equation for k 
In this section, we derive the equation for k. 

Let X : sp — » Quad be the Lie algebra isomorphism defined by Q(d,k,l) = 
I(S(d,k,l)) in theorem 3 in [1], where sp is the infinite dimensional Lie algebra 
containing 

'd k 



S(d,k,l) 



I -d 



DERIVATION 



7 



and k and I are symmetric, and Quad is the Lie algebra consisting of quadratics of 
the form 

Q{d,k,l) :=i (a x a*) 



1 

~ 2 

Recall that 
1 



Bjo! + a^ota; 1 f 

d{x,y) dxdy+- / k(x,y)a x a y dxdy 

l(x, y)a*a* dx dy . 



j!6L A, [A, [A, [A,V]]]] 



v{x -y,y-z) \4>{x)\ 2 {(j){y)(j>{z)a v a z + <j){y)(j){z)a* y a* z )dxdydz 
v(x — y,y — z) \<j){x)\ 2 <j)(y)4>{z)a* z a v dxdydz 



\ i 2 2 

+ 2 / v ( x -y>y- z )m x )\ my)\ a* z a z dxdydz 



we write 



n °t ) and M = ( L n 

—5 / V — TO 



with 



9 = -A5 X2 - (J v 12 ,23 I </> 3 1 2 dz^j^x^ Ul2 .23 I <^*2 1 2 1 03 1 2 dydz^jS 12 



and 

TO= ^ Wl2,23 I0 3 | 2d2: ^0102 

Though it's true that 

B=X(K) , 

for 

k(t,x,y) 



<"» "la,.) 

it's not quite true that 

ff + ii[A [A [A, [A,V]]]]=X(G + M) . 

However, the commutators of X(G + M) and ff + if [A, [A, [A, V}]}] with B 
are the same as in the discussion in page 15 in pQ. 
Now, Lq reads 

= H G +l( (l^-e K )e- K + [e K ,G]e- K + e K Me- K 



H G +I(Mi+M 2 + M 3 ) , 
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where 

Hq = Hq + y wi2,23 \4> 3 \ 2 <j>2(l>iala y dxdydz + ^ J v 12 ,23 |0 2 | 2 \<f> 3 \ 2 a*a x dxdydz, 
in this setting. 

Then by the definition of the isomorphism X, the coefficient of a x a y is —(.Mi, 12 + 
M.2,12 + ^3,12), and the coefficient of a*a* is (Mi, 21 + M.2,21 + -A^3,2i)- To write 
it explicitly: 

(4.2) - (Mi.u + M 2 ,u + M 3A2 ) 



= 0*1,21 +M 2 ,2l +^3,21) 
=(zsinh(/c) t + sinh(fc)g T + <?sinh(/c))cosh(/c)- 
(icosh(fc)i — [cosh(fc), <?]) smh(/s) 



— sinh(fc)msinh(fc) — cosh(fc)mcosh(fc) . 

Setting formula R~2"l to gives the equation l2.3l Hence we have the following theorem 
similar to Corollary 1 in [T]. 

Theorem 3. If (j> and k solve Wl^ and WIA then the coefficients of a x a v and a* a* 
drop out and Lq becomes 

Lq =H + J v 12 ,23 |0 3 | 2 ^lO-laydxdydz + i / t> 12;2 3 |0 2 | 2 \4>?\ 2 a* x a x dxdydz 

d(t,x,y) y dxdy , 

where d is given by formula\1.4\ and the full operator reads 



L=H + J wi2 :23 |03 1 2 (j) 2 4>ia x a y dxdydz + i / wi 2:23 |0 2 | 2 103 1 2 a* x a x dxdydz- 

d(t 7 x,y)a;a x dx+-±- 2 e B Ve- B + -l^e B [A,V]e- B + ^e B [A, [A,V]] e - B 

+ -^-re B [A, [A, [A, V]]]e~ B - N Xo ~ Xl 
36iV 2 

■ =L- N Xo - Xi , 

and 



Xo = ^ I v(x-y,y- z) \(j}(x)\ 2 \<j)(y)\ 2 \4>(z)\ 2 dxdydz, 
Xi(t) =-- I d(t,x,x)dx . 



This matches equation \2.4 
Notice that 



(^ B Ve- B + lN^e B lA lV }e- B 



-Ln-^A, [A,V]]e- B + ^N^e B [A, [A, [A, V]}} e - B jn , 

in our case. This ends the proof of theorem O 

The verification that the hypotheses of the main theorem are satisfied by some 
singular potentials will be provided in a subsequent article. 
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